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Abstract:  We  propose  a  novel  framework  for  a  new  class  of  two-channel  biorthogonal  filter  banks.  The 
framework  covers  two  useful  subclasses:  (i)  causal  stable  HR  filter  banks;  (ii)  linear  phase  FIR  filter  banks. 
There  exists  a  very  efficient  structurally  perfect  reconstruction  implementation  for  such  a  class.  Filter  banks 
of  high  frequency  selectivity  can  be  achieved  by  using  the  proposed  framework  with  low  complexity.  The 
properties  of  such  a  class  are  discussed  in  detail.  The  design  of  the  analysis/synthesis  systems  reduces  to  the 
design  of  a  single  transfer  function.  Very  simple  design  methods  are  given  both  for  FIR  and  IIR  cases.  Zeros  of 
arbitrary  multiplicity  at  aliasing  frequency  can  be  easily  imposed,  for  the  purpose  of  generating  wavelets  with 
regularity  property.  In  the  IIR  case,  two  new  classes  of  IIR  maximally  flat  filters  different  from  Butterworth 
filters  are  introduced.  The  filter  coefficients  are  given  in  closed  form.  The  wavelet  bases  corresponding  to  the 
biorthogonal  systems  are  generated.  We  also  provide  a  novel  mapping  of  the  proposed  one  dimensional  (ID) 
framework  into  two  dimensional  (2D).  The  mapping  preserves:  (i)  perfect  reconstruction;  (ii)  stability  in  the 
IIR  case;  (iii)  linear  phase  in  the  FIR  case;  (iv)  zeros  at  aliasing  frequency;  (v)  frequency  characteristic  of 
the  filters. 
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1.  INTRODUCTION 

Fig.  1.1(a)  shows  a  two-channel  maximally  decimated  filter  bank,  and  Fig.  1.1(b)  shows  the  well-known 
polyphase  form  for  this  system.  The  applications  of  such  multirate  systems  are  well-known  [1]— [7].  If  for 
all  input  x(n),  the  output  of  the  system  x(n)  =  cx(n  —  no)  for  some  nonzero  constant  c  and  integer  no,  the 
system  is  called  a  perfect  reconstruction  (PR)  system.  In  the  maximally  decimated  case,  PR  is  equivalent 
to  biorthogonality  [5].  A  number  of  PR  or  nearly  PR  systems  have  been  reported  before.  In  this  paper  we 
develop  several  new  results  for  two-channel  biorthogonal  filter  banks  based  on  a  useful  class  of  polyphase 
matrices. 

1.1.  Previous  work 

In  FIR  filter  banks,  all  the  four  filters  Ho,  Hi,  Fo,  and  Fi,  are  FIR  filters  while  in  the  case  of  HR  filter 
banks,  some  or  all  of  these  filters  are  HR  filters.  The  earliest  good  designs  for  the  HR  case  were  such  that 
the  analysis  bank  was  paraunitary  and  the  polyphase  components  of  Hq(z)  and  H\{z)  were  allpass  (see  [7] 
or  pp.  201  of  [2]).  Even  though  all  the  HR  filters  are  causal  stable,  the  reconstructed  signal  suffers  from 
phase  distortion.  HR  PR  filter  banks  typically  have  noncausal  stable  filters  or  causal  unstable  filters  [8],  [9], 
[10],  Recently  the  authors  in  [11]  proposed  a  HR  PR  technique  providing  causal  stable  solutions,  but  no 
satisfactory  design  method  was  given. 

In  earlier  design  of  two  dimensional  (2D)  filter  banks,  separable  filters  have  been  considered  because 
of  their  advantage  of  low  complexity.  However  nonseparable  filters  offer  more  freedom  in  the  design  and 
hence  in  general  will  give  better  performance.  Recently,  some  results  on  the  nonseparable  filter  banks  have 
emerged  [12]— [14],  However,  few  design  techniques  are  available  for  nonseparable  PR  filter  banks.  In  [12], 
a  design  method  based  on  space  domain  approach  is  given.  In  [13],  a  subclass  of  2D  paraunitary  systems 
(which  can  be  represented  as  a  cascade  of  ID  paraunitary  systems  of  degree  one)  is  considered.  However  in 
both  of  the  polyphase  approaches  above,  the  optimization  in  the  designs  involves  a  large  number  of  nonlinear 
constraints.  Thus  other  approaches,  such  as  ID  to  2D  mapping,  have  been  considered  [14]  — [19],  In  [14], 
even  though  PR  property  is  preserved  by  the  mapping,  the  frequency  responses  of  the  filters  will  change. 
In  [15]  and  [16],  a  mapping  of  ID  filter  banks  to  2D  filter  banks  is  given.  The  authors  apply  the  technique 
on  a  ID  two-channel  orthogonal  HR  system  to  achieve  a  2D  HR  filter  bank.  The  resulting  systems  have 
either  phase  distortion  or  stability  problem.  In  [17],  the  authors  employ  McClellan’s  transformation  on  the 
ID  maximally  fiat  FIR  halfband  filters  to  obtain  a  2D  biorthogonal  filter  bank.  However  because  of  the  lack 
of  factorization  theorems  in  the  2D  case,  one  of  the  lowpass  filters  is  constrained  to  have  all  its  zeros  at 
the  aliasing  frequency.  And  there  is  no  simple  way  to  ensure  the  frequency  selectivity  of  all  the  filters.  In 
[18],  the  authors  introduce  a  mapping  which  can  be  viewed  as  the  generalization  McClellan’s  transformation. 
2D  two-channel  PR  systems  with  good  frequency  selectivity  can  be  obtained  by  judiciously  designing  the 
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mapping.  However,  the  mapping  works  for  the  FIR  case  only  and  the  resulting  filters  usually  have  a  large 
number  of  coefficients. 

Some  of  the  results  in  this  paper  were  reported  in  the  earlier  conference  papers  [20],  [21]  and  [22].  For 
the  ID  case  [20],  both  of  the  linear  phase  FIR  and  causal  stable  IIR  solutions  for  PR  filter  banks  similar 
to  those  proposed  in  this  paper  were  given.  For  the  2D  quincunx  case  [21],  the  authors  constructed  a  ID 
to  2D  mapping  (which  is  the  same  as  the  mapping  given  in  Section  5  in  this  paper)  which  preserved  many 
of  the  desired  properties.  However  many  of  the  properties  of  the  ID  and  2D  biorthogonal  systems  were  not 
addressed  in  [20]  and  [21],  for  example  the  problem  of  imposition  of  zeros  at  the  aliasing  frequency  which  is 
important  for  the  purpose  of  generating  smooth  wavelet  basis  functions. 

1.2.  The  new  idea  and  its  merits 

In  this  paper,  we  constrain  the  polyphase  matrix  E(z)  such  that  [detE(z)]  is  a  delay.  Furthermore  we 
consider  E(z)  and  R(z)  to  be  either  (i)  both  causal  stable  IIR  or  (ii)  both  FIR.  In  each  case,  the  following 
properties  can  be  simultaneously  satisfied. 

1 .  Perfect  reconstruction  is  preserved  structurally  and  the  structural  complexity  is  very  low. 

2.  All  analysis  and  synthesis  filters  are  designed  by  controlling  a  single  transfer  function  /3(z)  [allpass  in 
the  IIR  case,  and  Type  2  (i.e.,  odd  order  symmetric  linear  phase  FIR)  in  the  FIR  case].  So  the  design 
procedure  is  very  simple.  It  is  very  easy  to  design  (3(z)  so  that  all  filters  have  good  responses  (lowpass 
or  highpass  as  the  case  may  be). 

3.  In  the  IIR  case,  all  the  analysis  and  synthesis  filters  are  causal  and  stable. 

4.  In  some  applications  such  as  image  coding,  the  linear  phase  property  of  the  analysis  and  synthesis  filters 
is  desired.  In  the  FIR  case,  the  filters  are  exact  linear-phase.  In  the  IIR  case,  we  can  force  the  phase 
response  of  the  filters  to  be  nearly  linear  in  the  passband,  as  we  shall  explain  and  demonstrate. 

5.  The  lowpass  analysis  filter  Ho(z)  can  be  forced  to  have  arbitrary  number  of  zeros  at  uj  =  ir.  Furthermore 
the  lowpass  synthesis  filter  Fq(z)  is  guaranteed  to  have  the  same  number  of  zeros  at  7r  as  Hq(z).  In 
both  of  the  IIR  and  FIR  cases,  we  give  closed  form  expression  for  the  filter  coefficients  that  provide 
maximum  number  of  zeros  at  7 r. 

A  new  class  of  biorthogonal  wavelet  basis  functions  can  be  generated  from  the  above  filter  bank.  The 
regularity  property  can  be  directly  controlled  by  imposing  multiple  zeros  at  7r  as  desired.  In  the  IIR  case, 
since  all  filters  are  causal  (in  addition  to  being  stable),  the  basis  functions  are  all  causal.  In  the  FIR  case,  the 
linear  phase  property  ensures  symmetry  of  the  wavelets,  while  at  the  same  time  providing  a  simple  control 
on  regularity  (because  the  number  of  zeros  at  tt  is  trivially  controlled). 

A  ID  to  2D  mapping 

Furthermore,  we  also  provide  a  novel  mapping  of  the  proposed  ID  filter  banks  into  the  2D  quincunx 
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case,  preserving  all  the  desirable  properties.  In  particular: 

1.  The  perfect  reconstruction  property  is  preserved. 

2.  In  the  HR  case  all  the  analysis  and  synthesis  filters  remain  causal  and  stable.  In  the  FIR  case  the  linear 
phase  property  is  preserved. 

3.  Even  though  the  filter  bank  is  nonseparable,  the  complexity  is  that  of  a  separable  filter  bank,  growing 
linearly  with  the  filter  order. 

4.  The  frequency  response  supports  for  the  filters  are  the  diamond  and  diamond- complement  as  desired 
for  the  quincunx  case  [15],  [2],  Moreover  the  filter  frequency  responses  are  ensured  to  be  good  simply 
by  designing  the  ID  filter  having  a  good  frequency  response.  Any  desired  specifications  can  be  met  by 
designing  a  ID  transfer  function  /3(z)  appropriately  as  we  shall  demonstrate. 

5.  If  the  ID  lowpass  filter  Ho{z)  has  k  zeros  at  7r,  then  the  resulting  2D  lowpass  filter  will  have  its  ith 
order  total  derivative  equal  to  zero  at  (n,  n),  for  i  =  0,  1,  . . k  —  1.  See  Sec.  5  for  details. 

We  also  provide  a  design  example  to  show  that  the  mapping  can  be  easily  applied  to  any  dilation  matrix 
(i.e,  decimation  matrix)  with  determinant  2. 

Relation  to  other  results  in  the  literature 

All  the  designs  proposed  in  this  paper  are  based  on  a  single  class  of  polyphase  matrices,  to  be  described 
in  Sec.  2.  However  some  of  the  filter  banks  reported  by  other  researchers  are  related  to  our  work.  In 
[23],  the  authors  derive  a  class  of  biorthogonal  linear  phase  FIR  filter  bank  which  turns  out  to  be  a  special 
case  of  our  two-channel  framework.  In  the  HR  maximally  flat  halfband  case,  our  solution  is  different  from 
the  traditional  HR  Butterworth  design  and  has  approximately  linear-phase  in  the  passband.  In  the  FIR 
maximally  flat  halfband  case,  the  solution  agrees  with  the  classical  FIR  maximally  flat  design  [24].  But 
our  construction  is  different  from  those  in  [25]  and  [6]  since  the  analysis  filters  are  factors  of  maximally  flat 
halfband  filters  in  [25]  and  [6]  while  our  analysis  filters  are  themselves  maximally  flat  halfband.  The  2D 
mapping  proposed  earlier  in  [15]  and  [16]  is  different  from  ours  because  it  is  known  that  the  earlier  mapping 
will  not  preserve  the  PR  property  in  general. 

1.3.  Outline  of  the  paper  and  notations 

Our  presentation  will  go  as  follows:  In  the  next  section,  we  will  derive  a  framework  for  the  two-channel 
biorthogonal  filter  banks.  Some  properties  of  such  class  will  be  described  in  detail.  In  Sec.  3,  we  will  discuss 
both  the  HR  and  FIR  filter  banks  which  are  covered  in  the  proposed  framework.  In  Sec.  4,  wavelet  basis 
functions  generated  from  the  proposed  filter  banks  will  be  presented  and  imposition  of  zeros  at  aliasing 
frequency  will  be  considered.  Two  new  classes  of  HR  maximally  flat  solution  are  given  in  closed  form.  In 
Sec.  5,  we  will  first  introduce  a  novel  2D  mapping  for  the  quincunx  case.  Some  properties  of  the  mapping 
are  discussed.  Then  both  the  HR  and  FIR  cases  are  considered.  Furthermore  numerical  examples  will  be 
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provided  throughout  the  discussion  to  demonstrate  the  idea. 

Notations  and  definitions:  Capital  boldfaced  letters  are  used  to  denote  matrices.  I  represents  the 
identity  matrix.  The  determinant  of  the  matrix  A  is  denoted  by  [det  A] .  Superscript  20  is  used  to  represent 
the  2D  function  obtained  by  applying  the  mapping,  for  example,  E2X>(z0,zi)  is  obtained  by  applying  the 
2D  mapping  to  E(z).  The  z-transform  of  h(n)  is  represented  by  H[z).  The  relation  between  the  filters 
{Hk(z),  Fk(z)}  and  the  polyphase  matrices  E(z)  and  R(z)  can  be  described  as  follows: 

Hk{z)  =  Ek,o(z2)  +  z-'Ek^z2),  and  Fk(z)  =  z_1f2o)fc(z2)  4-  R1}k(z2), 

where  Eij(z )  and  Rij(z)  are  respectively  the  ij-th  elements  of  the  matrices  E(z)  and  R(z).  A  filter  Hk{z) 
is  halfband  if  either  one  of  its  polyphase  components  Ekfi{z),  Ek,\(z)  is  a  delay. 


2.  A  FRAMEWORK  FOR  ID  BIORTHOGONAL  FILTER  BANKS 

Consider  Fig.  1.1,  where  a  two-channel  system  is  shown.  In  general,  R(z)  =  E-1(z)  for  perfect 
reconstruction.  It  is  not  easy  to  constrain  [detE(z)]  to  be  minimum  phase  for  stability  of  R(z),  so  let  us 
make  it  a  delay.  An  example  is 


-  (z~* 

V  0 


With  this  we  obtain 


H0(z)  =  z~2N  +  z~1/3(z2), 


but  H\{z)  =  z-(.2N'+t)  which  is  a  delay.  Thus  even  though  H0(z)  can  be  designed  to  be  a  good  lowpass  filter 
(as  we  will  show),  H\(z)  is  allpass  and  this  is  not  useful  for  subband  coding  applications.  We  can  modify 
ffi(z)  without  affecting  Hq(z)  by  taking  the  polyphase  matrix  to  be 


E (z)_(  0.5  0  \(z~N  f3(z)  \  _  (  0.5 z~N 

E(z)  ~  V  — 0.5a(z)  1 J  V  0  z-N  )-{  -0.5  z~Na 

Then  we  get  the  following  expressions  for  the  analysis  filters: 


0.5  (3{z) 

(z)  -0.5  a(z)f3(z)  +  z~N 


H0(z)  = 


(z~2N  +  z~l(3{z2)) 


,  Hx{z)  =  -a(z2)H0(z)  +  z~ 


Obtaining  ideal  responses  with  (2.4) 

First  notice  that  the  filter  Hq{z)  can  be  made  an  ideal  lowpass  filter  if  (3{z)  has  the  following  magnitude 


and  phase  responses: 


\P(ej2“)  1  =  1,  Vo, 


(2.5o) 


/  nr  j2u\  f  (-2 N  +  1  )w,  for  w  e  [0,7t/2]; 

^  '  \  (— 2AT  +  l)w  ±  7r,  for  w  €  (7r/2, 7r]. 


From  (2.4),  we  see  that  in  the  high  frequency  region,  Hi(eju)  has  unity  gain  since  \HQ{e^)\  =  0.  The 
function  a(z)  does  not  affect  Hq(z)  and  can  be  freely  chosen  to  shape  the  response  of  H\{z).  It  should 
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be  chosen  such  that  in  the  low  frequency  region,  a(z2)Ho(z )  cancels  with  z~2N  -1.  For  exact  magnitude 
cancellation,  |a(e,'“)|  must  be  unity.  Since  Hq(z)  is  linear  phase,  it  is  necessary  that  a(z)  has  linear  phase 
in  the  low  frequency  region.  Comparing  these  two  requirements  and  the  conditions  in  (2.5),  we  realize  that 
(3(z)  is  a  suitable  candidate  for  a(z).  Indeed,  if  N'  =  2 N  -  1,  H\{z)  is  an  ideal  highpass  filter.  In  this  case, 
we  have  an  ideal  filter  bank,  and  the  polyphase  matrix  E (z)  in  Fig.  1.1(b)  is 


E  (z)  = 


{  0.5  0  \fz~N  /?(*)  \_f 

\-o.5p(z)  iyV  o  z~2N+1 )  V-o 


0.5  z~N  0.5/3(z) 

).5z~n/3(z)  -0.5  (32{z)  +  Z~2N+1 


With  this  we  get  the  following  expressions  for  the  analysis  filters,  which  we  will  repeatedly  use  in  this  paper. 


H0(z) 


(. z~2N  +  z~l(3{z2)) 


,  H1(z)  =  -/3(z2)H0(z)  +  z-m+1. 


The  perfect  reconstruction  can  be  achieved  by  choosing  R(z)  in  Fig.  1.1(b)  to  be: 

RM_f*-2JV+1  -P(z)\  (  1  0  \_fz-2N+1-0.5p2(z)  -0.5 (3(z) 

1 \  0  Z~N  )  \0. 5 0(z)  0.5 )  V  0.5 z~Nf3(z)  0.5 z~N 

The  corresponding  synthesis  filters  can  be  verified  to  have  the  following  form: 


F0(z)  =  -ffi(-z),  F^z)  =  Ho(-z). 


This  choice  of  synthesis  filters  in  (2.9)  ensures  that  {Fo(z),  Fi(z)}  is  a  lowpass/highpass  pair  if  {Ho(z),  Hi(z)} 
is  a  lowpass/highpass  pair.  From  (2.6)  and  (2.8),  we  have  the  implementation  of  the  filter  bank  shown  in 
Fig.  2.1.  The  structure  is  similar  to  a  ladder  network  structure  [26], 

Remarks:  Of  course,  the  a(z)  in  (2.3)  can  be  taken  as  functions  different  from  P(z),  as  in  the  case  of 
[20],  [21],  [23].  This  will  provide  more  freedom  in  the  design.  However,  by  taking  them  to  be  the  same,  the 
biorthogonal  systems  can  have  some  additional  useful  properties.  Therefore,  we  will  only  consider  the  case 
when  a(z)  =  (3{z ). 

Two  useful  approximations  of  (2.5) 

The  ideal  choice  of  f3(z)  as  in  (2.5)  requires  infinite  complexity.  Therefore,  we  have  to  design  (3{z)  to 
approximate  the  conditions  in  (2.5).  However  the  approximation  will  not  change  the  perfect  reconstruction 
property  because  E(z)  in  (2.6)  and  R(z)  in  (2.8)  satisfy  R(z)E(z)  =  0.5z-3JV+1I,  regardless  of  the  choice 
of  f3(z).  Fig.  2.1  shows  that  the  frequency  responses  of  all  the  analysis  and  synthesis  filters  depend  on  one 
single  function  (5(z)  only.  The  frequency  selectivity  of  all  four  filters  depends  on  how  well  /3(z)  approximates 
conditions  (2.5).  This  makes  the  design  procedure  simple.  In  the  next  section,  we  will  provide  two  simple 
but  useful  approximations  which  correspond  to  the  following  two  cases: 

1.  Stable  HR  case:  Here  0(z)  is  chosen  to  be  a  causal  stable  allpass  function  so  that  (2.5a)  is  met  exactly. 
We  design  the  phase  response  of  the  allpass  filter  so  that  (2.5b)  is  approximately  satisfied.  This  leads 
to  a  biorthogonal  system  with  causal  stable  HR  analysis  and  synthesis  filters. 


5 


2.  Linear  phase  FIR  case:  To  satisfy  the  condition  (2.5b),  /3(z)  can  be  chosen  as  a  Type  2  linear  phase 
function  [2]  (filter  with  a  symmetric  impulse  response  of  even  length).  The  magnitude  response  of  /3 (z) 
is  optimized  to  be  as  close  to  unity  as  possible  so  that  (2.5a)  is  well-approximated.  This  leads  to  a  linear 
phase  biorthogonal  system. 

Additional  properties  of  the  filter  banks  designed  as  above 

In  Sec.  1,  we  have  outlined  some  properties.  Properties  1—4  mentioned  at  the  begining  of  Sec  1.2  are 
clear  from  the  above  discussion  and  Property  5  will  be  discussed  in  the  Sec.  4.  In  addition  to  these  five 
properties,  we  have 

1.  Double  halfband  property:  In  all  the  previous  constructions  of  two-channel  PR  filter  banks,  Ho(z)F<}(z) 
is  a  halfband  filter,  where  Ho(z )  is  not  necessary  halfband  but  a  factor  of  a  halfband  filter.  However  in 
our  construction  above,  one  can  verify  that  not  only  the  product  Hq(z)Fq(z)  but  also  the  filter  Ho{z) 
is  halfband. 

2.  Poles  of  filters:  In  the  HR  case,  notice  from  Fig.  2.1  that  there  is  no  feedback  loop  in  both  the  analysis 
and  synthesis  ends  in  the  ladder  network.  Therefore  the  filters  have  the  same  poles  as  those  of  (3(z 2) 
and  stability  depends  solely  on  the  allpass  function  /?(z).  Moveover  in  the  HR  case  if  the  allpass  filter 
/ 3{z )  is  implemented  by  using  the  robust  lattice  structure  [2],  the  filter  bank  is  stable  even  when  it  is 
realized  with  finite  wordlength. 

3.  Robustness  to  round  off  noise:  The  ladder  structure  shown  in  Fig.  2.1  is  similar  to  the  structure 
considered  in  [26].  By  using  the  same  reasoning  in  [26],  it  can  be  verified  that  the  round  off  noise  in 
the  analysis  end  is  compensated  by  that  in  the  synthesis  end.  Combining  this  with  the  structurally  PR 
property,  we  conclude  that  the  implementation  in  Fig.  2.1  preserves  PR  even  when  all  the  coefficients 
are  quantized  to  a  finite  percision  and  all  the  intermediate  results  are  rounded  off.  However,  if  the 
subband  signals  are  quantized  (which  is  usually  the  case),  this  property  is  lost. 

4.  Zeros  of  the  filters:  We  can  verify  that  F0(z)  and  Hffz)  in  (2.9)  and  (2.7)  can  respectively  be  rewritten 
as: 

F0(z)  =  (2 z~2N+1  -  f3(z2))HQ(z),  Hffz)  =  (2 z~2JV+1  +  /3(z2)) Fffz) .  (2.10) 

These  factorizations  give  the  filter  bank  an  interesting  structure  shown  in  Fig.  2.2.  From  (2.10),  it  is 
clear  that  if  j3{z)  is  FIR,  the  zeros  of  H0{z)  are  also  zeros  of  F0(z).  Even  when  (3(z)  is  an  irreducible  IIR 
transfer  function,  this  is  true  since  Ho(z)  is  in  the  form  of  (2.7)  and  the  zeros  of  denominator  of  /3(z2) 
cannot  cancel  the  zeros  of  Hq{z).  Moreover  if  \j3{e^)\  <  2,  both  Fo(eJ“)  and  Ho(eJ“)  have  the  same  set 
of  zeros  on  the  unit  circle.  The  same  is  true  for  the  pair  of  H\(z)  and  Fi(z).  In  particular,  if  Ho(z)  has 
r  zeros  at  z  =  -1,  this  implies  that  Fo(z)  has  no  fewer  than  r  zeros  at  the  same  point.  This  property 
is  important  in  the  generation  of  wavelets  since  for  biorthogonal  wavelets,  we  need  both  of  the  analysis 
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and  the  synthesis  wavelets  to  be  regular.  By  increasing  the  number  of  zeros  of  Ho(z)  at  z  =  —  1,  our 
construction  ensures  that  Fq(z)  has  at  least  the  same  number  of  zeros  at  2  =  -1.  This  is  the  property 
which  does  not  appear  in  the  previously  existing  constructions  of  biorthogonal  filter  banks. 

5.  Ripple  sizes  of  the  filters:  Since  H0{z )  is  a  halfband  filter  and  Hq(z)  +  F\(z)  =  z~2N,  we  have  the 
following  relationship  between  the  passband  ripple  Sp  and  the  stopband  ripple  8S : 

=  fi.(ffo)  =  6P(Fi)  =  S^Fj).  (2.11) 

Moreover  by  using  (2.10)  and  the  fact  that  f3(z 2)  w  — z_2Ar+1  in  the  high  frequency  region,  we  get 

6s(F0)k36s(H0),  *,(ffi)  »  3«.(Fi),  (2.12) 

(20 log 3  «  9.5  dB).  This  property  ensures  that  by  designing  Hq(z)  to  have  sufficiently  high  stopband 
attenuation,  we  can  ensure  that  all  the  other  three  filters  will  also  have  good  frequency  selectivity. 

6.  Complexity:  From  Fig.  2.1,  it  is  very  clear  that  the  analysis  and  synthesis  banks  have  the  same 
complexity.  Assume  that  fi(z)  has  order  TV.  For  the  HR  case,  by  using  the  one  multiplier  lattice 
structure  for  allpass  function  [2],  we  need  approximately  2 TV  multiplications,  6 N  additions,  and  57V 
delays.  For  the  FIR  case,  by  exploiting  the  symmetry,  we  need  approximately  TV  multiplications,  2 TV 
additions  and  3.57V  delays.  All  the  operations  are  at  a  lower  rate.  So  the  analysis  (or  synthesis)  bank 
requires  TV  and  0.57V  multiplications  per  input  sample  for  the  HR  and  FIR  case  respectively. 

7.  Near  linear  phase  in  the  HR  cases:  From  (2.7),  since  in  the  passband  the  magnitude  response  of  Hq(z) 
is  approximately  one,  the  transfer  function  (3(z2)  «  z~2N+1.  Therefore  Ho(z)  has  approximately  linear 
phase  in  the  passband.  Similar  arguement  is  true  for  H\{z). 


3.  DESIGN  PROCEDURES  FOR  THE  TWO  CLASSES  OF  BIORTHOGONAL  FILTER  BANKS 

In  this  section,  we  will  discuss  the  two  cases  of  the  approximations  of  (2.5)  given  in  the  last  section. 
Simple  design  procedures  will  be  given  for  both  cases. 

3.1.  Causal  Stable  IIR  Biorthogonal  Filter  Banks 

In  this  section,  /3(z)  in  (2.6)— (2.9)  is  taken  to  be  the  causal  stable  real  allpass  function 


ANl  (2) 


EiVl  — k 

k= 0  aNi,Ni-kz 

EJVi  „—k  * 

fc=0  aNi,kZ 


(3.1) 


where  qjv1i0  =  1  and  ^  are  real.  In  this  case,  Hq(z)  is  a  sum  of  a  delay  and  an  allpass  function.  See  Eq. 
(2.7).  It  is  an  IIR  halfband  filter  and  has  been  studied  by  some  researchers  [27],  [28].  Hq(z)  can  be  made 
lowpass  with  large  stopband  attenuation  and  small  passband  ripples  by  designing  the  phase  response  of  the 
allpass  function  to  approximate  (2.5b)  [29]. 
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Choice  of  N\:  From  the  monotone  decreasing  phase  property  [2]  of  a  causal  stable  allpass  function,  we 
know  that  the  phase  of  An1(z2)  spans  a  range  of  4iVi7r  when  u>  spans  a  range  of  2ir.  But  from  (2.5b),  /3(z2) 
spans  a  range  of  4iV7r  or  4 (N  —  l)7r.  To  make  the  range  spanned  by  both  of  the  functions  equal,  we  set 
JVi  =  N  or  N  —  1  and  this  results  in  two  classes  of  causal  stable  HR  filter  banks.  Since  the  derivation  and 
properties  of  both  of  the  classes  are  very  similar,  in  the  rest  of  the  paper,  we  consider  only  the  case  N\  =  N 
(we  will  point  out  at  those  places  where  the  second  class  has  a  different  property).  With  this  choice,  the 
analysis  filters  can  be  written  as 

H0(z)={z~-N.+\lAN-iz2)\  Hx(z)  =  -An(z2)H0{z)  +  z-4N+1.  (3.2) 

The  relationship  between  the  synthesis  and  analysis  filters  is  the  same  as  (2.9). 

Additional  properties  of  the  above  HR  filter  banks: 

1.  Preservation  of  zero  at  aliasing  frequency:  Substituting  z  =  —  1  into  the  expression  of  Hq(z)  in  (3.2), 
we  find  that  Hq(z)  always  have  a  zero  at  z  =  —  1,  independent  of  the  coefficients  ajv,fc.  In  particular, 
the  zero  is  preserved  even  when  all  a^,k  are  quantized  coarsely.  This  means  that  one  zero  at  2  =  —  1 
is  structurally  imposed.  This  is  important  in  the  generation  of  wavelet  bases  since  one  zero  at  z  =  —  1 
is  a  necessary  condition  for  the  existence  of  the  wavelet  functions  [6],  Note  also  that  Hi(z)  will  always 
have  a  structurally  imposed  zero  at  z  =  1 . 

2.  Low  sensitivity.  Since  there  exists  low  sensitivity  lattice  structure  for  allpass  function  [2],  the  filters 
have  low  passband  sensitivity.  Since  the  halfband  property  of  H0(z)  is  structurally  imposed,  it  has  low 
stopband  sensitivity  as  well. 

3.  Bump  in  the  transition  band:  Substituting  u >  =  n/2  into  the  expression  for  Hi(e?u)  and  Fo(e^u)  and 
using  the  fact  that  Ajv(— 1)  =  (— 1)N,  we  find  that  \Hi(ejul)\  =  |Fo(eJU,)|  =  y/2£  at  u  =  7r/2,  indepen¬ 
dent  of  the  allpass  function  An(z).  This  means  that  |£Ti(e,u’)|  and  |Fo(e-J1J)|  always  have  a  bump  of 
approximately  4  dB  at  w  =  7r/2,  no  matter  how  we  design  A^{z).  The  width  but  not  amplitude  of  the 
bump  can  be  reduced  by  increasing  the  complexity  of  A^{z). 

Example  3.1  ID  causal  stable  HR  filter  banks:  In  this  example,  N  =  3.  So  A^{z)  is  a  third  order 
allpass  function.  The  filter  bank  has  very  low  complexity:  To  implement  the  analysis  (or  synthesis)  bank, 
we  need  only  3  multiplications  per  input  sample!  By  using  the  eigenfilter  approach  for  allpass  functions 
[29]  we  optimize  the  coefficients  an  such  that  maximum  attenuation  in  the  stopband  of  Hq(z)  is  achieved. 
The  coefficients  are  obtained  as  03,1  =  0.473,  03,2  =  —0.094,  and  03,3  =  0.025.  For  the  filter  Ho{z),  the 
passband  edge  wp  =  0.47T  and  the  stopband  egde  u>3  =  0.67T.  The  stopband  attenuation  Ss(Hq)  =  41.9  dB. 
The  magnitude  responses  of  the  all  four  filters  are  shown  in  Fig.  3.1(a).  From  the  plots,  relations  of  ripple 
sizes  in  (2.11)  and  (2.12)  can  be  verified  and  it  is  clear  that  Ho(z)  and  Fq(z)  have  the  same  set  of  zeros  on 
the  unit  circle.  The  bump  of  approximately  4  dB  around  7r/2  is  clearly  seen.  The  group  delay  for  Hq(z)  and 


8 


Hx(z)  is  shown  in  Fig.  3.1(b).  The  filters  are  approximately  linear  phase  in  the  passband  and  the  stopband. 


3.2.  Linear  Phase  FIR  Biorthogonal  Filter  Banks 


In  the  linear  phase  FIR  case,  since  Hq(z)  is  a  linear  phase  halfband  filter,  it  can  be  designed  by  employing 
the  trick  developed  in  [30],  viz,  by  taking  f3(z)  in  (2.5) -(2.8)  to  be  a  Type  2  filter  [2]  which  has  a  symmetric 
impulse  response  of  length  2Ni.  In  this  case,  the  number  of  multiplications  required  to  implement  fi(z)  is 
N\,  the  same  as  the  iVith  order  allpass  function  An1  (z)  in  (3.1).  More  precisely,  let  Q(z)  have  the  following 
form: 

Vi 

V(z)  =  J2VkX  (z~Nl+k  +  z~Nl~k+1),  (3.3) 


k=l 


where  the  coefficients  vk  satisfy: 


N! 


sjT/vk  =  0.5, 

fc=  l 


(3.4) 


so  that  V(ej0)  =  1  and  7/o(ej0)  =  1.  It  is  well-known  that  a  Type  2  linear  phase  filter  always  has  a  zero  at 
z  =  —1.  In  order  to  satisfy  the  condition  (2.5b)  exactly,  it  can  be  verified  that  N\  should  be  equal  to  N.  By 
employing  the  trick  in  [30],  the  coefficients  vk  can  be  optimized  such  that  the  amplitude  response  of  V (eJ1J) 
is  as  close  to  unity  as  possible.  In  this  case,  the  analysis  filters  are: 


H0(z) 


(z~2N  +  z~lV(z2)) 
2 


Hi(z)  =  -V(z2)H0(z)  +  z~4N+1. 


(3.5) 


Example  3.2.  ID  linear  phase  FIR  filter  banks:  N  =  6.  To  implement  the  analysis  bank,  we  need  6 
multiplications  per  input  sample,  double  the  number  in  Example  3.1.  The  Type  2  linear  phase  function  V(z) 
is  designed  by  using  McClellan-Park  algorithm.  The  coefficients  are  obtained  as  =  0.630,  V2  =  —0.193, 
u3  =  0.0972,  v4  =  —0.0526,  vs  —  0.0272  and  v6  =  —0.0144.  For  the  filter  Hq(z),  the  passband  edge  u>p  =  0.47T 
and  the  stopband  edge  uis  =  0.67 r,  same  condition  as  Example  3.1.  The  stopband  attenuation  Ss(Hq)  =  39.2 
dB  and  6s(Hi)  =  30  dB.  The  magnitude  responses  of  all  four  filters  are  shown  in  Fig.  3.2.  The  relations  of 
ripple  sizes  in  (2.11)  and  (2.12)  can  be  verified. 

For  comparison,  we  will  consider  Johnston’s  design  with  nearly  the  same  specifications.  The  Johnston’s 
filter  24C  in  Appendix  7.1  of  [31]  has  6S  =  30  dB  and  ws  =  0.5867T.  For  Johnston’s  filter  32D,  8S  =  38  dB  and 
w,  =  0.5867T.  To  implement  the  analysis  bank,  we  need  respectively  12  multiplications  and  16  multiplications 
per  input  sample  for  the  above  two  cases.  Thus  as  compared  to  6  multiplications  in  our  filter  bank,  the 
Johnston’s  design  has  more  complexity  than  our  design.  Moreover,  there  is  reconstruction  error  (0.1  dB  for 
24C  and  0.025  dB  for  32D)  in  Johnston’s  filter  bank. 


4.  IMPOSITION  OF  MULTIPLE  ZEROS  AT  7T 

The  relation  between  continuous- time  wavelet  and  discrete-time  perfect  reconstruction  filter  bank  is  well 
known.  A  way  to  construct  the  scaling  and  wavelet  functions  from  the  filter  coefficients  was  first  given  by 
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Daubechies  in  [6].  Starting  from  the  impulse  response  coefficients  ho(n)  and  h\(n),  a  pair  of  continuous-time 

functions  cf>H0  (x)  and  (x)  are  constructed  such  that  they  satisfy: 

oo 

<pH0(x)  =  ^2hQ{n)<t>Ho{2x-n),  (4.1a) 

n=0 

oo 

Y^h1(n)cj)Ho(2x  -n).  (4.16) 

n=0 

Here  <j>Ha(x)  an<i  are  respectively  called  the  analysis  scaling  and  wavelet  functions.  For  the  synthesis 

end,  we  can  write  similar  expressions  for  the  synthesis  scaling  and  wavelet  functions,  <j>F0(x)  and  ipF1  ( x ).  The 
conditions  for  the  existence  of  such  limit  functions  were  given  in  [6],  It  is  always  desirable  to  have  smooth  or 
“regular”  limit  functions.  It  was  shown  that  in  order  to  achieve  limit  functions  of  high  regularity,  we  need 
to  have  a  sufficient  number  of  zeros  at  the  aliasing  frequency  tt.  Therefore  in  the  rest  of  this  section,  we  will 
show  how  to  impose  zeros  at  tt  for  the  proposed  filter  banks. 

4.1.  Causal  Stable  IIR  Wavelet  Bases 

For  the  purpose  of  achieving  regularity,  we  impose  multiple  zeros  of  Hq(z)  at  tt.  Since  the  denominator 
does  not  provide  any  zeros,  we  consider  only  the  numerator  of  Hq(z).  Except  for  a  delay,  the  numerator  of 

Ho{z)  can  be  written  in  terms  of  aw>fc  as  follows: 

N 

Pr(u)  =  ^ajv.fc  cos(2fc  -  l/2)w.  (4.2) 

k=0 

To  obtain  r  zeros  at  z  =  —1,  we  set 

P«)(7r)=^)Pfl(w)L=.  =  °>  for  *  =  1,  2,  •  •  • ,  r  -  1.  (4.3) 

Note  that  when  i  is  even,  p£'  (tt)  is  always  equal  to  zero.  This  proves  that  Pr(u)  always  has  an  odd 
number  of  zeros  at  w  =  it.  Therefore,  we  can  write  r  =  2ro  +  1.  In  this  case,  we  obtain  a  set  of  ro  linear 
constraints  as  follows: 

N 

“  4fc)2l_1  =0,  for  i  =  1,  2,  . . . ,  r0.  (4.4) 

fc=o 

The  set  of  linear  constraints  in  (4.4)  can  be  satisfied  exactly  in  the  optimization  of  the  phase  response  of  the 
allpass  function  An(z)  by  using  the  efficient  eigenfilter  approach  [29],  [32]. 

Maximally  flat  IIR  wavelets 

To  obtain  a  maximally  fiat  solution,  i.e.,  maximum  possible  number  of  zeros  at  tt  consistent  with  the 
constraint  that  Hq(z)  =  0.5 (z~2N  +z~1AN(z2))}  we  set  r0  in  (4.4)  as  large  as  possible.  However  if  r0  >  JV+1, 
then  we  can  list  the  first  (N  +  1)  linear  constraints  given  by  (4.4)  as  follows: 


V  andermonde 
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where  xk  =  1  -  4 k.  Since  all  the  xfc  are  nonzero  and  distinct,  the  two  matrices  on  the  left  hand  side  are 
nonsingular  and  hence  invertible.  We  get  [ajv,o  aN,i  ■  ■  ■  =  0  which  violates  the  requirement  that 

ajv,o  =  1.  This  proves  under  the  constraint  that  Hq(z)  =  0.5[z~2Ar  +  z~1An(z2)\,  the  filter  Hq(z)  can  have 
at  most  2 N  +  1  zeros  at  it.  Indeed  we  can  show  that  the  maximally  flat  HR  filter  has  exactly  2N  +  1  zeros 
at  7 r.  To  see  this,  we  set  tq  =  N  and  rewrite  the  set  of  N  linear  equations  given  by  (4.4)  as  follows: 


where  =  1  -  Ak  and  the  fact  that  aN0  =  1  has  been  imposed.  These  equations  fully  determine  A^(z) 
(hence  all  the  filters)  and  there  is  no  further  parameter  to  be  optimized  numerically.  As  the  matrices  are 
invertible,  the  solution  for  ajv.fc  always  exists  and  it  is  unique.  Furthermore,  it  is  shown  in  Appendix  A  that 
ajv,fc  has  the  following  closed  form  solution: 


O/V.fc 


(-l)fc-i  /N\  "  (2 «  -  1) 

2k-  1  \k)  (2A  +  2t-  1)’ 


0  <  k  <  N, 


(4.7) 


where  (^)  —  (N^ly~k\  ■  The  frequency  responses  of  Hq(z)  corresponding  to  N  =  1,  3,  •  •  • ,  10  are  shown  in 
Fig.  4.1.  Note  that  although  these  filters  have  a  numerator  of  degree  4AT  -  1  (excluding  the  trivial  delay 
factor),  they  have  only  2 N  +  1  zeros  at  z  ~  —  1.  This  implies  that  some  of  the  zeros  are  not  at  z  =  -1 
for  N  >  1  and  therefore  these  HR  maximally  flat  filters  are  different  from  the  Butterworth  halfband  filters. 
Moreover  they  have  nearly  linear  phase  in  the  passband,  as  justified  at  the  end  of  Sec.  2  and  demonstrated 
in  Fig.  3.1(b).  For  the  case  of  N  =  1,  one  can  verify  that  the  solution  is  a  third  order  Butterworth  filter. 


Remarks: 


1.  If  the  function  (3(z)  is  taken  as  (N  —  l)th  order  allpass  filter  (i.e.  Ni  =  N  —  1),  then  we  will  get  a 
second  class  of  causal  stable  HR  wavelet.  In  this  case,  under  the  constraint  that  Ho{z )  =  0.5[z-2Ar  + 
z_1Aiv-i(^2)],  the  process  of  imposition  of  zeros  at  7 r  is  very  similar  to  the  derivation  above.  The 
maximally  flat  HR  filter  Hq(z)  of  this  second  class  will  have  2 N  —  1  zeros  at  7r.  The  closed  form  solution 
for  ajv-i,fc  is  given  as: 


(_l)fe-l  /N  -  !\  T-r1  (2»  +  l) 

2k +  1  \  k  J  y  (2Jfe  +  2t  +  l)’ 


1  <k<N-l, 


(4.8) 


and  o  =  1. 


2.  Notice  that  for  a  perfect  reconstruction  system,  if  we  interchange  the  analysis  and  synthesis  filters,  the 
perfect  reconstruction  property  is  retained.  In  many  applications  such  as  coding,  compression,  storage 
and  approximation,  the  regularity  of  the  synthesis  functions  is  more  important  [17].  Thus  we  can  choose 
the  wavelet  with  higher  regularity  among  ipHo(x)  and  ipFo(x)  as  the  synthesis  wavelet. 
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Example  4-1-  We  generate  the  limit  functions,  <pHa>  i,  4>F0y  and  ipp1 ,  corresponding  to  the  filter  bank 
in  Fig.  1.1(a).  To  generate  the  analysis /synthesis  scaling  and  wavelet  functions,  we  use  the  cascade  algorithm 
in  [6]  for  eight  iterations.  We  consider  the  following  two  cases 

(i)  No  linear  constraint  is  set,  Ho(z)  has  only  one  zero  at  ir.  The  analysis  and  synthesis  filters  are  the 
same  as  those  in  Example  3.1.  For  the  analysis  bank,  the  scaling  and  wavelet  functions,  4>h0  and  i pglt 
are  respectively  shown  in  Fig.  4.2(a)  and  (b).  The  scaling  and  wavelet  functions  corresponding  to  the 
synthesis  bank,  (j>F0  and  V’Fi,  are  shown  in  Fig.  4.2(c)  and  (d). 

(ii)  As  a  comparison,  we  also  generate  the  scaling  and  wavelet  functions  corresponding  to  the  HR  maximally 
flat  filters  (^mai  and  ipmax)  for  N  =  3.  In  this  case,  the  filter  Hq(z)  has  seven  zeros  at  7r.  The  limit 
functions  are  shown  in  Fig.  4.3.  For  a  better  comparison  on  smoothness,  in  Fig.  4.4  we  show  a  zoom-in 
for  Fig.  4.2(a)  and  Fig.  4.3(a).  We  see  that  the  limit  functions  in  Fig.  4.3  are  more  regular  than  the 
functions  shown  in  Fig.  4.2. 


4.2.  Linear  phase  FIR  wavelet  bases 

To  impose  multiple  zeros  at  it  for  the  linear  phase  FIR  case,  the  procedure  is  very  similar  to  that  given 
above.  Another  set  of  linear  constraints  can  be  obtained  and  incorporated  in  the  procedure  of  optimization. 
It  can  be  verified  that  for  this  case,  Hq(z)  always  has  an  even  number  of  zeros  at  it. 

Maximally  flat  linear  phase  FIR  wavelets 

The  FIR  maximally  flat  filters  have  been  studied  by  a  number  of  researchers  [24],  [33],  [6],  [25],  In 
[6]  and  [25],  a  maximally  flat  halfband  FIR  filter  is  used  to  construct  compactly  supported  maximally  flat 
wavelets.  In  our  linear  phase  FIR  filter  bank,  if  all  the  freedom  is  used  to  impose  zeros  at  it,  we  will  arrive 


at  the  same  solution  as  that  in  [6]  and  [25].  The  closed  form  solution  for  FIR  maximally  flat  halfband  filters 
was  in  [33],  [25]  as: 

_  (-pw^natw+i/z-i)  . 

or\r  u\\r\T  i  i  uwtou  i\‘  ^  ' 


K  2(N  —  k)\(N  —  1  +  k)\(2k  —  1) '  v  ' 

Differences  between  our  construction  and  those  in  [25],  [6]:  In  [25],  Hq(z )  is  taken  to  be  a  factor 
of  a  maximally  flat  halfband  filter.  In  [6],  power  spectral  factorization  is  considered.  However,  in  our  linear 
phase  structure,  Hq(z)  is  taken  to  be  this  halfband  filter  itself,  and  not  a  factor.  Since  the  H'o(z)  constructed 
in  [6]  is  a  power  spectral  factor  of  the  Ho(z)  in  our  structure,  our  linear  phase  scaling  function  <j>Lp{x)  is 
related  to  that  constructed  by  Daubechies  in  [6],  <Pd{x)  as: 


4>Lp{x)  =  <j>D{x)  *  <P*d{-x)> 


(4.10) 


where  *  denotes  convolution  and  <j>*D  denotes  the  complex  conjugate  of  q>D-  From  (4.10),  it  is  clear  that  the 
regularity  of  4>lp(x)  is  twice  that  of  4>d(x).  However  the  order  (and  the  number  of  zeros  at  it)  of  Hq(z)  in 
our  construction  is  twice  that  of  Hq(z)  in  the  construction  in  [6],  Comparing  the  complexity,  both  of  the 
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constructions  have  approximately  the  same  number  of  multiplications  (because  in  our  construction,  linear 
phase  property  can  be  exploited) . 

Example  J^.2.  In  this  example,  we  construct  the  limit  functions  corresponding  to  the  filter  bank  in  Fig. 
1.1(a)  for  the  linear  phase  FIR  case.  The  cascade  algorithm  is  used  for  eight  iterations.  We  consider  two 
cases 

(i)  First,  Hq(z)  is  designed  such  that  no  linear  constraint  other  than  (3.4)  is  satisfied,  therefore  it  has  two 
zeros  at  n.  The  analysis  and  synthesis  filters  are  the  same  as  those  in  Example  3.2.  The  limit  functions 
(4>h0i  V’ff'D  $Fq  and  i/'fJ  are  respectively  shown  in  Fig.  4.5(a),  (b),  (c)  and  (d). 

(ii)  For  a  comparison,  we  show  the  limit  functions  of  the  maximum  flat  case  (<£max  and  ipmax)  for  N  =  6. 
In  this  case,  Hq(z)  has  twelve  zeros  at  n.  The  plots  are  shown  in  Fig.  4.6.  It  can  be  verified  that  the 
limit  functions  in  Fig.  4.6  are  smoother  than  those  in  Fig.  4.5. 

5.  MAPPING  INTO  2D  QUINCUNX  PERFECT  RECONSTRUCTION  FILTER  BANKS 

In  this  section,  we  will  generalize  the  ID  framework  discussed  in  Sec.  2  to  the  2D  case.  We  will 
focus  on  the  quincunx  subsampling  case  which  has  the  subsampling  lattice  shown  in  Fig.  5.1.  Notice  that 
the  dilation  matrix  has  determinant  2.  The  corresponding  maximally  decimated  filter  bank  has  only  two 
channels.  Furthermore  it  represents  the  simplest  nonseparable  subsampling  lattice. 

In  the  2D  case,  we  know  that  the  desired  passband  supports  of  the  filters  depend  not  only  on  the  lattice 
but  also  on  the  choice  of  dilation  matrix  M  [12].  In  the  rest  of  this  section,  we  will  consider 


The  coset  vectors  are  respectively: 

k°=(s)'  k>=0)'  (5-2> 

With  this  M,  the  ideal  supports  for  alias  free  decimation,  SPD(nM~T)  [Chapter  12,  2]  is  shown  in  Fig. 
5.2,  where  the  diamond  and  diamond-complement,  fio  and  fix,  correspond  to  the  low  frequency  and  high 
frequency  regions  respectively.  One  can  verify  that  M  defined  in  (5.1)  has  its  eigenvalues  A*  equal  to  ±\/2 
and  M2  =  21.  It  has  a  dilation  in  both  the  directions.  Therefore,  M  satisfies  the  conditions  for  a  well-behaved 
matrix  defined  in  [19].  Given  the  dilation  matrix  M  as  in  (5.1)  and  the  coset  vectors  in  (5.2),  the  simple 
delay  chain  system  and  the  noble  identities  are  shown  in  Fig.  5.3(a)  and  (b)  respectively. 

Although  the  discussion  in  this  paper  is  mainly  on  the  quincunx  subsampling  case  with  the  dilation 
matrix  M  and  the  coset  vectors  ki  defined  above,  we  will  provide  a  design  example  in  the  last  section  to 
show  that  the  method  discussed  in  this  section  can  be  easily  generalized  to  any  2D  system  with  decimation 
matrix  M  having  [detM]  =  2. 
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5.1.  A  ID  to  2D  Mapping 


In  this  subsection,  we  will  first  give  a  2D  mapping  and  then  apply  the  mapping  to  the  framework 
developed  in  Sec.  2.  Given  any  ID  biorthogonal  systems  with  the  polyphase  matrices  of  the  form  in  (2.6) 
and  (2.8),  we  will  use  the  following  transformation  on  the  polyphase  components: 

(i)  First  replace  the  ID  transfer  function  /3(z)  with  the  separable  2D  transfer  function  /?(zo)/3(zi). 

(ii)  Replace  all  the  remaining  ID  delay  z~1  with  the  2D  delay  zj^zf1. 

This  results  in  nonseparable  analysis  and  synthesis  filters  as  we  will  see.  Under  this  transformation,  the 
polyphase  matrices  E2D(zo,  z{)  and  R2D(z0,  zi)  of  the  2D  system  can  be  written  respectively  as: 

0.5  0\  / {z0z1)-N  (3(zo)fi(zi) 

-0.5/3(zo)/3(zi)  1/  \  0  (z0zi)~27V+1 

0.5(z02i)~JV  0.5/3(z0)f3(zi) 

„  -0.5(z0zi)~n/3(z0)/3(-zi)  -0.5/32(z0)/?2(zi)  +  (z0zi)-2JV+1 , 


E2U(zo,z1)  = 


(5.3) 


'(20*1) 


-2JV+1 


R (zo,zi)  = 


~fi(zo)P{zl)\ 


(5.4) 


0  (zozi)-N  )  \0 .5/3(zp)/3(zi)  0.5y 

r  (^02i)”2Ar+1  -  0.5/32(z0)/?2(z1)  — 0.5/3(zo)/3(zi)  ' 

O.5(z0zi)-'iv/3(zo)/3(zi)  0.5(z0zi)“n 

From  the  above  two  equations,  we  have  the  implementation  of  the  2D  perfect  reconstruction  filter  bank  as 
Fig.  5.4.  By  using  the  noble  identities  in  Fig.  5.3,  we  can  write  the  analysis  and  synthesis  filters  as: 

Z02N  + 


Hq{z0,zi)  = 

#1(20,  *i)  =  -f3{zozT1)Pizoz\)Ho{zo,zi)  +  zo4N+1, 


(5.5a) 


Fb(zo,  zi)  =  —Hi(—zo,  -zr),  Fi(z0,zi)  =  H0{-z0,-zi).  (5.56) 

Comparison  of  the  above  transformation  with  those  in  [15]  —  [17]:  McClellan’s  transformation  is 
used  in  [17]  to  obtain  a  FIR  maximally  flat  halfband  filter.  The  transformation  proposed  in  this  paper  differs 
from  McClellan’s  transformation  in  the  sense  that  the  former  operates  on  the  polyphase  components  while 
the  latter  operates  directly  on  the  filter.  In  [15]  and  [16],  the  authors  obtain  a  2D  filter  bank  from  ID  by 
employing  the  following  transformation: 

E?f(z0,z1)  =  EiJ(zo)EiiJ(z1),  (5.6) 

where  E ij  is  the  (i,  j)th  element  of  E.  It  is  clear  that  in  our  transformation  E^(z0,  zi)  ^  Eiii(zo)Ejii(zi). 
Therefore  our  mapping  is  different  from  that  in  (5.6). 
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Properties  of  the  proposed  2D  mapping 

Properties  1—5  in  Sec.  2  continue  to  hold  after  minor  modifications  to  suit  the  2D  context.  In  addition, 
the  2D  filter  bank  satisfies  the  following  properties 

1.  Double  halfband  property:  It  is  easy  to  see  that  Hq(zq,  z\)  satisfies  #o(zo,  z\)  +  Hq{—zq,  —zi)  =  zf2N 
and  Ho(zo, zi)Fo(zo, zi)  satisfies  a  similar  property.  This  is  the  extention  of  the  ID  double  halfband 
property  in  the  2D  quincunx  case. 

2.  Stability  of  the  2D  analysis  and  synthesis  filters:  If  the  ID  transfer  function  /3(z)  is  causal  then  so 
are  the  functions  f3(zo)f3(zi)  in  Fig.  5.4.  That  is  /3(zo)/5(zi)  is  a  first-quadrant  filter  (the  impulse 
response  is  zero  unless  no  >  0  and  ni  >  0).  If  (3{z)  is  BIBO  stable,  then  so  is  fi(zo)j3(zi)  so  that  the 
polyphase  matrix  in  Fig.  5.4  is  also  BIBO  stable.  Since  the  analysis  filters  are  obtained  from  this  stable 
structure,  these  filters  are  guaranteed  to  be  BIBO  stable.  However  we  see  that  the  term  /3(zozf *)  has 
entered  the  expressions  for  the  analysis  filters  because  of  the  noble  identities,  see  Fig.  5.3(b).  It  can  be 
shown  that  this  violates  the  condition  for  the  so-called  first-quadrant  stability  (p.  166  of  [35]).  This  is 
explained  by  the  fact  that  the  analysis  filters  are  not  first-quadrant  filters,  even  though  BIBO  stable. 
This  is  consistent  with  the  observation  that  the  quincunx  decimator  M  in  (5.1)  has  the  negative  entry 
-1.  Indeed,  the  expression  y(n)  —  i(Mn)  means  y{n0,ni)  =  x{n0  +  n\ , n0  —  rij)  so  that  there  is  a 
time-reversal  operation  buried  in  the  decimation  process.  The  same  remarks  apply  for  the  synthesis 
filters,  that  is  the  2D  synthesis  filters  Fq{zo,zi)  and  F\(z0,zi)  are  BIBO  stable  even  though  they  are 
not  first-quadrant  filters. 

3.  Perfect  reconstruction  is  preserved. 

4.  If  the  ID  lowpass  filter  Hq(z)  has  k  zeros  at  7r,  then  the  frequency  response  of  7/o(eJW°,  e^1 )  can  be 
written  as 


Ho (eJW° ,  eJWl )  =  (l  +  e^^)fcP1(c;o,Wi)  -  (l  +  Wl),  (5.7) 


where  |jPi(vt,  tt)  |  and  \P'2(r,  w)  |  are  finite  quantities.  The  proof  of  (5.7)  is  given  in  Appendix  B.  Notice 
that  both  of  the  factors  [1  +  e-°-5J(wo+“i)]  an(j  jq  e-o.5j(wo-wi+2w)j  are  zero  at  (7r,7r).  Furthermore 
one  can  verify  that  all  the  mixed  partial  derivatives  satisfy 

tori+'<L  (58) 

From  (5.8),  we  conclude  that  the  total  derivatives  [36] 


dnH0{n,ir)  =  J2  Qdwj'K 


dn 


doJodu’f  1 


—  0,  for  n  <  k. 


(5.9) 


According  to  [34],  (5.9)  is  a  necessary  condition  for  the  regularity  of  2D  wavelet.  The  necessary  and 
sufficient  condition  is  still  unknown. 
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5.  In  the  FIR  case  the  linear  phase  property  of  the  analysis  and  synthesis  filters  is  preserved. 

6.  In  the  HR  case,  the  2D  analysis  and  synthesis  filters  have  a  line  of  zeros  in  the  frequency  plane  at  wo  =  0 
or  at  wo  =  7T. 


Proof:  Substituting  zq  =  —1  into  the  expression  for  Hq(zq,  z\)  in  (5.5a)  and  using  the  fact  that  P{zqz\) 
is  allpass,  one  immediately  finds  that  Ho(-l,zi)  =  0,  V  z\.  Since  Fb(zo,zi)  contains  Hq(zq,  z\)  as  a 
factor,  Fo(-l,  z\)  =  0.  Similarly,  we  can  prove  that  Fi(l,  z\)  =  H\(l,  zi)  =  0,  V  z\. 

7.  The  lowpass/highpass  characteristics  of  the  frequency  responses  of  the  filters  are  preserved. 

Proof:  Assume  that  /3(z)  satisfies  the  ideal  conditions  in  (2.5).  Then  we  have 


/?(z0zi)  = 


(z0zi)  2  ,  for  €  [0,  7t/2]; 

~(zozi)~^d±,  for  €  (tt/2, vr]. 


(5.10a) 


P(z0z1  *)  = 


1  f  (z0z11) 


2 


for  **=&>  €  [0, 7r/2] ; 


(5.106) 


\  -(zo^i  for  g  (vr/2, 7r]. 

By  using  the  above  equations,  we  find  that  ^(zozi)^(zozj-1)  is  equal  to  zq2N+1  when  (wo,  wi)  6  Ho  and 
equal  to  —  Zq2N+1  when  (wo,  wj)  e  Hj.  This  proves  that  Ho(zo,  zi)  has  the  ideal  diamond  support  Ho- 
Similarly  it  can  be  shown  that  Hx(zo,  zi)  will  have  the  support  of  ideal  diamond-complement.  Thus 
when  the  conditions  in  (2.5)  are  well-approximated  by  the  ID  transfer  function  /3(z),  the  response  of 
the  2D  filters  will  be  good. 

Comments  on  the  complexity:  Though  the  2D  analysis  and  synthesis  filters  are  nonseparable,  it  is  clear 
from  the  expressions  for  the  polyphase  matrices  that  the  complexity  is  comparable  to  that  of  a  separable 
filter  bank.  More  precisely,  it  is  equal  to  twice  the  complexity  of  the  ID  transfer  function  /3(z). 


5.2.  2D  Nonseparable  Filter  Banks 

In  this  section,  we  will  give  two  numerical  examples  to  demonstrate  the  mapping  proposed  above.  We 
separately  apply  the  2D  mapping  to  the  filter  banks  in  Example  3.1  (IIR)  and  Example  3.2  (FIR)  respectively. 

Example  5.1.  2D  IIR  filter  banks:  In  this  example,  we  transform  the  ID  filter  bank  in  Example  3.1  into 
the  2D  case  by  using  above  mapping.  Since  N  =  3,  the  allpass  function  As(z)  needs  only  3  multiplications. 
Since  the  complexity  of  the  2D  analysis  (or  synthesis)  bank  is  equal  to  twice  that  of  As(z),  we  need  only  6 
multiplications  per  input  pixel  to  implement  the  analysis  (synthesis)  bank.  The  responses  of  Ho(zo,zi)  and 
Hi(*b,*0  are  shown  in  Fig.  5.5(a)  and  (b)  respectively.  The  supports  of  the  two  filters  are  diamond  and 
diamond-complement  respectively  as  desired.  The  stopband  attenuation  5s(Ho)  ~  42  dB  and  6s(Hi)  «  32 
dB.  Again,  we  see  that  Hi  is  about  10  dB  worse  than  Ho  in  the  stopband.  The  line  of  zero  of  Hi  at  wo  =  0 
is  clearly  seen  in  Fig.  5.5(b). 


Example  5.2.  2D  FIR  filter  banks:  In  this  example,  the  ID  filter  bank  in  Example  3.2  is  transformed 
into  the  2D  case.  To  implement  the  2D  analysis  (or  synthesis)  bank,  we  need  12  multiplications  per  input 


pixel.  The  magnitude  responses  of  Ho(zo,z\ )  and  H\(zo,zi)  are  shown  in  Fig.  5.6(a)  and  (b)  respectively. 
The  stopband  attenuation  6s(Hq)  «  40  dB  and  S3(H\)  «  30  dB. 

6.  CONCLUDING  REMARKS 

In  this  paper,  we  have  derived  a  framework  for  a  new  class  of  two-channel  biorthogonal  filter  banks. 
The  filter  banks  under  the  framework  allow  a  structurally  perfect  reconstruction  implementation  as  in  Fig. 
2.1.  It  is  interesting  that  we  can  arrive  at  precisely  the  same  ladder  in  Fig.  2.1  by  using  the  novel  approach 
in  [26]  developed  for  a  totally  different  application,  namely  cancellation  of  roundoff  error.  The  proposed 
systems  have  very  low  complexity.  Filter  banks  of  high  frequency  selectivity  can  be  achieved  by  controlling 
a  single  transfer  function  (3{z)  in  Fig.  2.1.  Two  different  choices  of  /3(z)  lead  to  causal  stable  HR  and  linear 
phase  FIR  filter  banks  respectively.  The  properties  of  the  proposed  filter  banks  were  discussed  in  detail.  We 
showed  that  zeros  at  aliasing  frequency  can  be  imposed.  Two  new  types  of  HR  maximally  flat  filters  were 
derived  and  the  solutions  were  given  in  closed  form.  In  addition  to  perfect  reconstruction  property,  these  HR 
filters  have  nearly  linear  phase  in  the  passband.  Furthermore,  we  also  mapped  the  ID  filter  banks  derived  in 
this  paper  into  2D  cases.  The  design  of  a  2D  biorthogonal  (stable  HR  or  linear  phase  FIR)  filter  bank  reduces 
to  the  design  of  a  single  ID  transfer  function.  The  new  transformation  preserves  many  of  the  properties  of 
the  ID  systems.  Before  we  conclude  the  paper,  we  would  like  to  provide  an  example  to  demonstrate  that 
the  mapping  in  Sec.  5  can  be  easily  generalized  to  arbitrary  dilation  matrix  M  with  determinant  equal  to  2. 

Example  6.1.  2D  HR  filter  banks:  The  ID  prototype  filter  bank  is  taken  to  be  that  in  Example  3.1.  The 
dilation  matrix  and  the  coset  vectors  are  respectively: 

M=  (l  O'  ko=(o)  “>-©•  <6'1) 

With  the  above  matrix  and  coset  vectors,  the  ideal  passband  support  for  the  Ho{zq,z\)  is  SPD(-nM.~T), 
which  is  shown  in  Fig.  6.1  (shaded  area).  By  using  the  tranformation  introduced  in  Sec.  5,  we  find  that 
the  polyphase  matrices  in  this  example  are  the  same  as  those  in  Example  5.1.  Thus  it  also  has  very  low 
complexity.  The  only  differences  are  the  dilation  matrix  and  the  coset  vectors.  With  the  M  and  ki  chosen 
as  (6.1),  the  responses  of  Hq(zq,zi)  and  H\{zq,  z{)  are  shown  in  Fig.  6.1(a)  and  (b).  We  see  that  Ho  and 
H i  have  approximately  the  desired  support. 


17 


Appendix  A.  Proof  of  Eq.  (4.7) 
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Fig.  1.1.  (a)  A  two-channel  analysis/synthesis  filter  bank,  (b)  A  redrawing  of  (a)  by  using  the  polyphase  repre¬ 
sentation. 

Fig.  2.1.  An  implementation  of  the  proposed  biorthogonal  filter  bank. 

Fig.  2.2.  A  redrawing  of  Fig.  2.1,  where  Hq(z)  =  0.5 {z~2N  +  z-1/J(z2)),  and  Fj(z)  =  Ho(-z). 
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analysis  and  synthesis  filters,  (b)  Group  delays  of  Hq(z)  and  H\{z). 

Fig.  3.2.  (Example  3.2)  Magnitude  responses  of  the  linear  phase  FIR  filter  bank. 

Fig.  4.1.  Magnitude  responses  of  the  IIR  maximally  flat  filters  of  the  form  0.5 [z~2N  +  z_1Aat(«2)]  where  An{z) 
is  a  iVth  order  allpass  function,  for  N  =  1,  2,  . . . ,  10. 

Fig.  4.2.  (Example  4. l(i))  The  limit  functions  generated  by  using  the  IIR  filter  bank  in  Example  3.1  (f/o(z) 
has  one  zero  at  x).  (a)  Analysis  scaling  function,  (b)  Analysis  wavelet  function,  (c)  Synthesis  scaling 
function,  (d)  Synthesis  wavelet  function. 

Fig.  4.3.  (Example  4.1(h))  The  limit  functions  generated  by  using  the  IIR  maximally  flat  filter  bank  (HQ(z)  has  7 
zeros  at  it),  (a)  Analysis  scaling  function,  (b)  Analysis  wavelet  function,  (c)  Synthesis  scaling  function, 
(d)  Synthesis  wavelet  function. 

Fig.  4.4.  A  zoom-in  for  Fig.  4.2(a)  and  Fig.  4.3(a)  demonstrating  the  improved  “regularity”  obtained  by  imposing 
zeros  at  n. 

Fig.  4.5.  (Example  4.2(i))  The  symmetric  limit  functions  generated  by  using  the  FIR  filter  bank  in  Example  3.2 
(H0(z)  has  2  zeros  at  t r).  (a)  Analysis  scaling  function,  (b)  Analysis  wavelet  function,  (c)  Synthesis 
scaling  function,  (d)  Synthesis  wavelet  function. 

Fig.  4.6.  (Example  4.2(h))  The  symmetric  limit  functions  generated  by  using  the  FIR  maximally  flat  filter  bank 
(H0(z)  has  12  zeros  at  tt).  (a)  Analysis  scaling  function,  (b)  Analysis  wavelet  function,  (c)  Synthesis 
scaling  function,  (d)  Synthesis  wavelet  function. 

Fig.  5.1.  The  quincunx  subsampling  lattice. 

Fig.  5.2.  Ideal  supports  for  alias-free  decimation  in  quincunx  case. 

Fig.  5.3.  Some  details  for  the  quincunx  decimator.  (a)  The  delay  chain  perfect  reconstruction  system,  and  (b) 
the  noble  identities. 

Fig.  5.4.  A  2D  biorthogonal  filter  bank  obtained  from  Fig.  2.1  by  mapping. 

Fig.  5.5.  (Example  5.1)  Magnitude  responses  of  the  perfect  reconstruction  IIR  analysis  bank,  (a)  Ho(z0,  z\)  and 
(b)  Hi(z0,  zi).  The  normalized  frequency  =  Vi/in. 

Fig.  5.6.  (Example  5.2)  Magnitude  responses  of  the  perfect  reconstruction  FIR  analysis  bank,  (a)  Hq{zq,  z\)  and 
(b)Fi(zo,zi). 

Fig.  6.1.  The  ideal  supports  for  alias-free  decimation  for  M  defined  in  (6.1). 

Fig.  6.2.  (Example  6.1)  Magnitude  responses  of  the  perfect  reconstruction  analysis  bank  with  the  decimator  M 
defined  in  (6.1).  (a)  Ho(zq,zi)  and  (b)  H\(zq,  z\). 
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Fig.  5.3.  Some  details  for  the  quincunx  decimator 

(a)  The  delay  chain  perfect  reconstruction  system,  and 

(b)  the  noble  identities. 
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Fig.  5.4.  A  2D  biorthogonal  filter  bank  obtained  from  Fig.  2.1  by  mapping. 
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Fig.  6.2.  (Example  6.1)  Magnitude  responses  of  the  perfect  reconstruction  analysis 
bank  with  the  decimator  M  defined  in  (6.1). 

(a)  Ho(zo,zi)  and  (b)  Fh(zo,zi). 
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